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ABSTRACT 

About half of the baryons in the local Universe are thought to reside in the so-called warm-hot 
intergalactic medium (WHIM) at temperatures of 0.1-10 million K. Thermal soft excess emission in 
the spectrum of some cluster outskirts that contains OVII and/or OVIII emission lines is regarded as 
evidence of the WHIM, although the origin of the lines is controversial due to strong Galactic and solar 
system foreground emission. We observed the Coma- 11 field, where the most prominent thermal soft 
excess has ever been reported, with Suzaku XIS in order to make clear the origin of the excess. We did 
not confirm OVII or OVIII excess emission. The OVII and OVIII intensity in Coma-11 is more than 
5(7 below that reported before and we obtained 2a upper limits of 2.8 and 2.9 photons cm~ 2 s~ 4 sr^ 1 
for O VII and O VIII, respectively. The intensities are consistent with those in another field (Coma-7) 
that we measured, and with other measurements in the Coma outskirts (Coma-7 and X Com fields 
with XMM- Newton). We did not confirm the spatial variation within Coma outskirts. The strong 
oxygen emission lines previously reported are likely due to solar wind charge exchange. 
Subject headings: clusters: individual: 
universe 
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1. INTRODUCTION 



Based on N-body simulatio ns of cosmological lar ge- 
scale structure formation fe.g. JCen fc Ostrikerll2006f) . a 
significant (30-50%) fraction of baryons are thought to 
reside in a 'warm-hot' gaseous phase (T = 10 5-7 K), 
which is hard to detect with the instruments currently 
in operation. This warm-hot gas, whose density is 10 -6 - 
10~ 4 cm -3 , is called the warm- hot intergalactic medium 
(WHIM) . Firm detection of the WHIM is important be- 
cause it is the most promising candidate for the "missing 
baryons"; i.e., it is thought to explain the discrepancy 
betwe en the baryon densit y observed in the local uni- 
verse ( Fukugita et al,lll998D and that in the distant uni- 
verse ( Rauch et al.lll997h or that calculated from the ob- 
served fluctuations o f the cosmic microwave background 
(jSpergel et alj 120071) . Although - 10% of the baryons 
has been resolved as the WHIM in the temperature range 
T = 10 5 ~ 6 K through O VI absorption features ob - 
served with FUSE and HST (jDanforth & Shulll 120071 ) . 
more baryons are thought to reside in a hotter phase at 
T = 10 6 ~ 7 K. The WHIM of this temperature may be 
detected via emission or absorption lines from highly ion- 
ized elements in X-ray spectra, such as O VII and O VIII. 

Despite many observations so far, the existence of the 
WHIM in the hotter phase is not yet confirmed (see re- 
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view of iBregmanl [2007. and references therein). Stud- 
ies of the WHIM using absorption line s only lead to 
one d isputed detection toward Mkn 421 (jNicastro et alJ 
12005ft . a result w hich was not confirmed later and was 
thus questioned (jRasmussen et al.l l2007t iKaastra et al.1 
|2006| ). Possible detections of emission from the WHIM 
have been reported in cluster outskirts, as a soft X-ray 
(E 5, 0-5 keV) excess containing O VII and O VIII lines 
IKaastra et al.ll2003tlFinoguenov et alll2003h . The origin 
of O line emission is, however, controversial due to strong 
Galactic and solar system foreground e mission, which has 
a spectral shape similar to the WHIM (|Bonamente et all 
l20MlBre gman &: Llovd-Dav ics 2006). In particular, X- 
rays induced by solar wind charge exchange (SWCX) 
could produce serious contamination because of i ts hig h 
variability (e.g., ICravensI 120001 : iFuiimoto et~aT1 120071 ). 
Note that there has been no detection of rcdshiftcd emis- 
sion lines from the WHIM so far. 

The most prominent thermal soft excess ever reported 
is in the Coma outskirts , especially in the Coma-11 field 
(jFinoguenov et al.ll2003l . hereafter FBH03). FBH03 also 
observed spatial variation within the Coma outskirts; 
Coma-7 shows roughly five times smaller intensity than 
Coma-11. Not only O lines, but also an EUV excess 
(|Bonamente et ai]|2003|) and ma rginal detection of Ne IX 
emission and absorption lines (jTakei et al.l I2007L here- 
after T07) are reported in the Coma outskirts. These 
results may support the existence of the WHIM. 

We observed the Coma-11 an d Coma-7 fields 
with the XIS l|Koyama et al.l l2007| ) onboard Suzaku 
(jMitsuda et al.ll2007f ). in order to make clear the origin of 
the excess. The XIS is an ideal instrument to study the 
emission lines from the WHIM, because it does not show 
a large low-energy tail in the pulse-height distribution 
function, even in the very soft (E < 0.5 keV) band. 

We assume a Hubble constant of 70 km s _1 Mpc -1 
or hjQ = 1 and fi m = 0.3, = -7. The so- 

lar metal abundance is given by lLoddersI ((2003). Un- 
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TABLE 1 
Observation Information 



TABLE 2 

Best-fit parameters of ComaBKG. 





Coordinates 


Exposure a 




(J2000.0) 




Coma- 11 


(12 h 58 m 31.7 s , 28°24'13.2") 


53.02 ks 


Coma-7 


(12 h 57 m 22.9 s , 28°08'58.1") 


25.26 ks 


ComaBKG 


(13 h 15 m 00.0 s , 31°39'28.0") 


30.78 ks 



Parameter 



Unit 



Value 



a After screening as described in the text 

less otherwise stated, errors and upper limits are at 
the 68 % confidence level in the figures, and at the 
90 % confidence level in the text and tables. The emis- 
sion line intensities are quoted in line units (LU) de- 
fined as photons cm~ 2 s _1 sr _1 . One LU corresponds 
to 8.46 x 10 -8 photons cm -2 s _1 arcmin~ 2 . 

2. OBSERVATIONS AND DATA REDUCTION 

We observed three fields from 2007 Jun 19 to 2007 
Jun 22; two of them are in the Coma outskirts, namely 
Coma- 11 and Coma-7 that FBH03 observed with XMM- 
Newton, and the other is an offset observation (hereafter 
we call it ComaBKG) of 5.0-degree distance from the 
Coma center (NGC 4874). The coordinates and net ex- 
posure times are summarized in Table [T] The XIS was 
operated in the normal mode and with the spaced-row 
charge injection (SCI). We analysed the clean-event data 
of Ver 2.0 processing with standard parameters, but with 
the latest gain calibration (Ver 2.6.1.16). We excised the 
last 6.7 ks of the Coma-7 observation, during which the 
XIS suffered a problem with zero-level determination. 

We excluded the regions where point sources are 
found or which arc irradiated by the Fe calibration 
sources. The point sources were identified using the 
XMM-Newton Serendipitous Source Catalogue (2XMM; 
IWatson et aTll2008f) . as well as by examining the XIS im- 
ages by eyes. Response files were created using Suzaku 
ftools xisrmfgen version 2007-05-14 and xissimarfgen ver- 
sion 2007-09-22 with CALDB released on 2008-01-14. We 
assumed uniform emission from the sky. The spectra and 
response files of the two XIS-FIs (XIS0 and 3) arc added 
using ftools mathpha, marfrmf, and addrmf. X-ray spec- 
tra are analysed with XSPEC 12.4.0. 

3. SPECTRAL ANALYSIS 
3.1. ComaBKG field 

ComaBKG is observed in order to precisely determine 
the foreground and background emission, in particular 
O and Ne line intensities around the Coma cluster. The 
instrumental background, estimated using a ftool xis- 
nxbgen ver sion 2007-11-23 w ith COR2 as the sorting 
parameter (|Tawa et al.ll2008T) . were subtracted from the 
ComaBKG spectra. 

We modeled the spectra with a sum of two thin ther- 
mal plasma models (APEC in XSPEC) and a power-law 
to represent emission from the local hot bubble (LHB), 
Milky- Way halo (MWH) and extragalactic cosmic X- 
ray background (CXB), respectively. MWH and CXB 
were convolved with absorption in the Galaxy (WABS 
in XSPEC). We fitted the spectra in the energy range 
0.5-5.0 keV for the XIS-FI and 0.35-5.0 keV for XIS1 
(XIS-BI). The spectra of XIS-FI and XIS-BI were fit- 
ted simultaneously. We ignored 1.40-1.55 keV where a 



Galactic Nh cm -2 
LHB temperature keV 
LHB emission measure pc cm" 

MWH temperature keV 
MWH emission measure pc cm" 

CXB photon index 
CXB normalization 3 
X 2 /dof thermal model 



O Vll intensity 
O VIII intensity 
Ne IX intensity 
Ne X intensity 
X 2 /dof with lines 



1.0 X 10 20 (fixed) 

0.07 (fixed) 
(3.0 ±3.0) X 10~ 2 

u - i44 -0.027 

(8-5±3.o) x 10-3 
1.52 ±0.09 
7.7 ±0.5 
153.40/113=1.36 



LU 
LU 
LU 
LU 



6.5 ± 1.0 
1 1+0.9 

15+ 015 
U - iO -0.14 

(< 0.2) 
140.22/105=1.34 



In a unit of photons cm" " 2 s" 1 keV" 1 sr" 1 at 1 keV 



0.01 




1 2 

Energy (keV) 

Fig. 1. — Spectra and best-fit models (top) and residuals in units 
of sigma (bottom) of ComaBKG. Red and black show BI and FI, 
respectively. Instrumental background is subtracted. 



relativel y large uncertaint y exists in the instrument back- 
ground (|Tawa et al.ll2008f ) The redshift and abundance of 
the two thermal plasma models were fixed to and 1.0 
solar, respectively. The temperature of the LHB was also 
fixed to 0.07 keV, while the temperature of the MWH was 
allowed to vary. The column density of t he Galactic hy- 
droge n gas was fixed to 1.0 x 10 20 cm -2 (jKalberla et alJ 
120051 ). Figure [Hand Tabic [2] show the best-fit model. 

We also fitted the spectra with another model in order 
to determine O VII, O VIII, Nc IX and Ne X intensity 
more precisely. Four Gaussians were added to represent 
the four line emissions, while O and Ne abundance of 
the LHB and MWH were fixed to zero. The energies of 
the Gaussians were allowed to vary ±35 eV around the 
resonance line energy (0.574 keV for O VII, 0.653 keV for 
O VIII, 0.922 keV for Ne IX and 1.022 keV for Ne X) in 
order to compensate the uncertainty in the energy-scale 
calibration, while we fixed the line width to zero. The 
intensities of the four lines arc also shown in Table [H 
The relatively large y 2 value may be due to incorrect 
calibration of the energy resolution. The \ 2 improved to 
128.18 when we added 50 eV (FWHM) additional width 
into the model with lines. The best-fit parameters were 
not sensitive to this calibration uncertainty; they stay 
the same within the error at 90% confidence level. We 
confirmed that the fitting is not improved when we allow 
N-r to vary. 



Lack of O-line excess in the Coma outskirts 



3 



TABLE 3 

Best-fit parameters of ComaII and Coma7. 



Parameter 


Unit 


Coma- 11 


Coma-7 


Galactic Nn 

Temperature 

Abundance 

Redshift 

EM 

X 2 /dof without lines 


_2 

cm 

keV 

solar 

pc cm~ 6 


1.0 X 10 (hxed) 

6 7 +0A 
°- '-0.3 

0.36 ±0.11 
0.0231 (fixed) 
0.100 ±0.002 
307.51/275=1.12 


1.0 X 10 (fixed) 
o O+0.8 
°- d -0.6 

0.21 ±0.15 
0.0231 (fixed) 
0.110 ±0.002 
295.40/281=1.05 


O Vll intensity 
O VIII intensity 
Ne IX intensity 
Ne X intensity 
X 2 /dof with lines 


LU 
LTJ 
LU 
LU 


0.0 (< 0.9) 
1.1 ±0.8 

0.1 (< 0.4) 
0.3 ±0.2 
297.71/271=1.10 


0.2 (< 2.1) 

1.7 ± 1.1 
0.4 (< 0.8) 
0.3 (< 0.6) 
284.92/277=1.03 




Energy (keV) Energy (keV) 

Fig. 2. — Spectra of Coma-11 (left) and Coma-7 (right). Red and black data correspond to BI and FI respectively. ComaBKG spectra 
were subtracted to correct the background emission. 



3.2. Coma-11 and Coma-7 fields 

The back- and foreground emission of Coma-11 and 
Coma-7 were corrected by subtracting the ComaBKG 
spectra that were extracted from the same region as 
Coma-11 or Coma-7 in the detector coordinates. Wc 
then fitted the spectra 0.5-7.5 keV for XIS-FI and 0.35- 
7.5 keV for XIS-BI) by the emission of the intraclus- 
ter medium (ICM), i.e., a thin thermal plasma model 
(APEC) convolved with Galactic absorption (wabs). The 
spectra of XIS-FI and XIS-BI were again fitted together. 
The spectra and the best-fit models are shown in Fig- 
ure [21 and the best-fit parameters are in Table [31 We did 
not see a soft excess, i.e., emission from the WHIM. Next, 
we added four Gaussians to the model, corresponding to 
lines of O VII, O VIII, Ne IX and Ne X at the cluster 
redshift, in order to investigate possible contributions of 
the WHIM. The central energy of each line was fixed to 
that of the resonance lines at redshift z = 0.0231, and the 
width of the line was fixed to zero. The results are shown 
in Tabled The fit is slightly improved (A% 2 ~ 10) by 
adding the lines. However, the significance of each line is 
less than 3cr. The fitting without FI also gave the same 
results within the statistical errors, which excludes the 
possible systematic bias due to the low statistics of the 
FI data below 1 keV. We also tried to fit the spectra 
with two thin thermal plasma models, i.e., one for the 
intracluster medium and one for the WHIM. This did 
not improve the fit. To obtain robust upper limits of line 
intensities, we consider systematic uncertainties in cali- 
bration: the centroid energy (±35 eV), the energy res- 



olution (10-80 cV FWHM), energy scale (±50 eV), and 
the contamination thickness (±5 x 10~ 17 cm" 2 ). The 2a 
upper limits we obtained for O VII, O VIII, Ne IX, and 
Ne X are 1.4, 2.7, 0.5, and 0.9 LU for Coma-11 and 2.6, 
3.7, 1.2, and 1.1 LU for Coma-7, respectively. The uncer- 
tainty in Galactic O intensity should also be taken into 
account. From the variation among background fields in 
T07, we adopted ±1.2 LU for O VII and ±0.6 LU for 
O VIII as la. Then, O VII and O VIII 2a upper limits 
are 2.8 and 2.9 LU for Coma-11 and 3.5 and 3.9 LU for 
Coma-7, respectively. 

4. DISCUSSION 

We measured the spectra of Coma-11 and Coma-7 with 
the Suzaku XIS. The spectra were fitted with a single 
thin thermal plasma model and the emission from the 
WHIM was not confirmed. Here we compare our results 
with previous XMM-Newton observations. We do not 
discuss the lower energy excess (E < 0.5 keV) observed 
with ROSAT and EUVE, because the XIS is not sensitive 
in the energy. 

FBH03 observed Coma-11 and Coma-7 with XMM- 
Newton. The O VII and O VIII intensity measured 
in FBH03 are 6.4 ± 0.8 and 6.2 ± 0.8 LU in Coma-11, 
and 1.4 ± 0.8 and 0.9 ± 0.6 LU in Coma-7, respectively. 
The intensities of Coma-7 in FBH03 are consistent with 
our Suzaku observation, while those of Coma-11 are dis- 
crepant at the 5a level. The best-fit ICM temperature 
and abundance in Coma-11 are also inconsistent. Even if 
we fixed the temperature and abundance to the best-fit 
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Fig. 3. — Solar wind proton flux as a function of time, around Coma-11 and Coma-7 observations with XMM-Newton (left) and around 
our observations with Suzaku (right), measured with ACE SWEPAM instrument. The time of XMM-Newton and Suzaku observations is 
indicated with vertical lines. 



values of FBH03, the discrepancies of the O line inten- 
sities remain as large. The O intensities observed with 
Suzaku are similar to those of Coma-7; we did not con- 
firm the spatial variation. 

T07 measured O VII, O VIII and Ne IX emission line 
intensities with XMM-Newton in the field of X Com and 
four background fields whose distance from NGC 4874 is 
46.2-162.7 arcmin. The intensities of O VII, O VIII and 
Ne IX lines measured by T07 are 6.6 ± 0.2, 1.8 ± 0.2, and 
0.7 ±0.1 LU in the X Com field and 7.4 ± 1.2, 2.6 ± 0.6, 
and 0.4 ±0.1 LU, in the (averaged) background fields, 
respectively. The values of the X Com field shown above 
contain contributions of the background emission. The 
intensities of the X Com field are consistent within 90% 
confidence range with those of the Coma-11 and Coma-7 
measured in this work, and those of the background field 
of T07 are also consistent with those of ComaBKG. 

The origin of the large discrepancy in Coma-11 O in- 
tensities between this work and FBH03 could be con- 
tamination by SWCX emission. The solar wind pro- 
ton flux observed with ACE (Advanced Composition Ex- 
plorer ) showed flares both during Coma-11 observation 
of FBH03 (2.5xl0 9 cm' 2 s" 1 ) and during ComaBKG 
observation of this work (1.5xl0 9 cm~ 2 s _1 ) as shown 
in Figure [31 Both arc ~ 8 times larger than the qui- 
escent state around each observation time. Hence, two 
possibilities should be considered: O intensities of Coma- 
11 in FBH03 are overestimated, or those of ComaBKG in 
this work are overestimated and thus those of Coma-11 in 
this work are underestimated. However, the latter case 
is not likely because the O-line intensities are consistent 
with or even smaller than T07 measurements of back- 
ground fields around the Coma cluster. There is no indi- 
cation of the geocoronal SWCX emission; the first point 
in the line of sight (LOS) where the geomagnetic field 
starte d open was always higher than 10 times Earth ra - 
dius 8 (jFuiimoto et al.ll2007tlTsvganenko fc Sitnovll2005D . 
and there was no variation in the O VII and O VIII 
flux during the ComaBKG observation. The heliosphcric 
SWCX emission could be weak since it is faint at so- 
lar minimum on the high ecliptic lattitude line of sight 
(jKoutroumpa et aTll2007f ). Moreover, even in the case 
that the ComaBKG observation was heavily contami- 

7 |http://www.srl.caltech.edu/ACE/ASCl 

8 Calculated using the software GEOPACK- 

2005 and T96 magnetic field model 

(http://geo.phys . spbu.ru/~tsyganenko/modeling.html 

), with CUAWeb s olar-wind parameters 

(http: / /cdaweb . gsfc .nasa.gov/ cdaweb/sp_phys/ ) 



nated with the SWCX emission, the 2a upper limits of 
Coma-11 O VII and O VIII intensity would be no larger 
than 4.8 LU and 3.5 LU, respectively, which are still 
lower than values of FBH03. These values are calculated 
by assuming the lower limits of O VII and O VIII inten- 
sities as the LHB contribution without any MWH con- 
tribution, measured toward a m olecular cloud MB M 12 
(3.5 LU and LU, respectively; ISmith et al.ll2007t ). On 
the other hand, it is more likely that the FBH03 observa- 
tion suffered heliosphcric SWCX contamination, since a 
similar solar wind flare during XMM-Newton observation 
of the Hubble deep field north increased the O VII and 
O VIII flux by 6 25±0.67 and 5.53±2.87 LU, re spectively 
(jSnowden et al.ll2004l ; iKoutroumpa et al.ll2007D . The so- 
lar cycle at the observations, the position of the Earth 
and the LOS, which also affects the X-ray flux induced by 
SWCX, are similar between the two observations. Note 
that the observation of FBH03 was performed at solar 
maximum, while our observation was at solar minimum. 
The field-to-field variation of O VII and O VIII flux re- 
ported by FBH03 might also be due to the variation of 
the SWCX flux. 

YT and SK thank T. Hagihara for his instruction on 
the calculation of the geomagnetic field distribution. YT 
is financially supported by the Dutch Organization for 
Scientific Research (NWO). EDM and JNB gratefully ac- 
knowledge support for Suzaku activities from NASA. 
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